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SUMMARY 


A computerized algorithm to generate cross-sectional dimensions and fiber 
orientations for composite airframe structures is described, and its application 
in a wing structural synthesis is established. The algorithm unifies computa- 
tions of aeroelastic loads, stresses, and deflections, as well as optimal struc- 
tural sizing and fiber orientations in an open-ended system of integrated computer 
programs. A finite-element analysis and a mathematical-optimization technique 
are the main components of the procedure. Design constraints include stress, 
strain, local buckling, and minimum gage. The algorithm is applied to evaluate 
a series of basic design alternatives such as types of construction, types of 
material, and manufacturing restrictions for the low-aspect-ratio wing of a large 
transport airplane. Included also is a review of the computational efficiency 
of the method . 


INTRODUCTION 

Extensive use of computer-aided design in synthesis of airframes for a low- 
aspect-ratio supersonic cruise airplane is reported in references 1, 2, and 3. 

The studies in reference 1 focused on the arrow-wing configuration. The objec- 
tive was to size the wing structure for minimum mass, subject to static strength 
requirements. The basic wing was assumed to be constructed of ribs, spars, and 
sandwich cover panels. The variants considered in the sizing studies were an 
all-metal (titanium) wing, a wing of hybrid construction in which the cover panels 
were made of a graphite-polyimide composite, and, finally, an all-composite wing. 

Because the airplane (ref. 1) was large and very flexible, it was important 
to include the effect of structural deformations on aeroelastic loads. This situ- 
ation required detailed aerodynamic and structural analyses which were computa- 
tionally very large. In some of the studies presented in reference 1, composite 
materials were assumed to be used for wing construction. Composite materials 
allow the designer the added flexibility of adjusting the proportions of mate- 
rial at the various orientations and the orientation angles themselves. This 
added design flexibility increases the number of design variables and further 
adds to the complexity of the sizing problem. 


As indicated in reference 4, most available procedures for optimizing 
composite structures are for isolated components of large structures. (See 
refs. 5 to 7.) In addition, there is a continuing effort for improvement in 
the efficiency of optimization procedures to extend their applicability to 
entire structures. (See refs. 8 to 13.) Because of computational size and 
need to optimize at the detail design level, the subject problem entails char- 
acteristics of the synthesis on both component and large structure levels. 

The purpose of this paper is to document the procedure used in the studies 
of reference 1, to explain the algorithm details, and to give a more detailed 
account of several structural sizing trade-offs. 



Nick J. Santoro, Southern Company Services, Inc., Birmingham, Alabama, 
participated in the conceptual part of the reported work while with Vought Cor 
poration, Hampton Technical Center, Hampton, Virginia, and was responsible for 
generation of the numerical results. 


SYMBOLS 

Values are given in both SI and U.S. Customary Units. The calculations 
were made in U.S. Customary Units. 

A,B,D matrices of membrane stiffness, bending plate stiffness, and bending 
membrane coupling stiffness, respectively 

a,b rectangular panel dimensions (see fig. 2) 

gi dimensionless measure of constraint satisfaction 

I set of constraints 

k vector of plate curvatures due to bending 

M vector of moments 

m mass of core face-sheet bonding material per unit area 

n»dep'"'indep nonoptimum structural mass dependent on panel area and indepen- 
dent of panel area, respectively 

mp structural mass of panel defined by equation (10) 

N vector of membrane force 

Nj{fNy,N}jy normal and shear membrane forces acting on edges of panel (see 

fig. 2) 

panel area 

thickness of material whose fibers are oriented at a given angle 
"caliper" thickness of sandwich panel (face-to-face) 
sum of thicknesses of both face sheets 

thickness of material whose fibers are oriented at (j) and y, 
respectively 

thickness of ith ply and total number of plies in one face sheet of 
symmetric sandwich 

behavior variable 


S 

t 
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tt 
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Ya limit on behavior variable (allowable value) 

Z defined by equation (6) 

e vector of membrane strains Gy, Gjjy or e-] , G 2 , g-|2 

p mass density of ply material 

Pj, mass density of core 

Ot stress in fiber 

P 2 stress in matrix, normal to Oi, in ply plane 

Pit allowable tension limit of fiber 

^2t'^2c allowable tension and compression matrix limits 
T12 in-plane shear stress 

(j),Y orientation angles of fibers (see fig. 2) 

Subscripts; 

a allowable limit 

cr critical 

Configurations; 

B metal baseline construction 

IM/E low Young's modulus/hybrid construction 

HM/H high Young's modulus/hybrid construction 

low Young's modulus/pure composite construction 
HM/P high Young's modulus/pure composite construction 

LM/H-CA conservative alldwable strain variant of LM/H configuration 

OVERALL DESIGN APPROACH 

Airplane flexibility introduces a coupling of aerodynamic loads and struc- 
tural displacements. This coupling can be accounted for by organizing the anal- 
yses involved into an iterative process. The flow chart in figure 1 (from 
ref. 14) depicts a particular organization of the iterative process used in 
this study. 
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Two iterative loops are shown in the flow chart of figure 1 . Loop I begins 
with the computation of aerodynamic loads on the entire airplane whose shape is 
preset to the defined cruise shape (best for cruise aerodynamic efficiency) . 
Initially, the same shape is used to compute loads for both cruise and maneuver. 
The loads are subsequently applied to the airframe using the initial cross- 
sectional dimensions, and the ensuing analysis yields displacement and stresses 
for both cruise and maneuver . The cruise displacement outputs are used to cal- 
culate the jig shape (shape to which the airframe is to be built to attain the 
desired shape under cruise loads) . The jig shape and the maneuver deflections 
are superimposed to calculate maneuver shape. The loop is closed by feeding back 
the maneuver shape to calculate a new set of aerodynamic loads for maneuver . 

In the flow chart of figure ^ , loop II involves use of the stress outputs 
to resize the airframe structural components. Resizing of structural components 
produces new cross-sectional dimensions (e.g., wing-cover thicknesses) which are 
fed back to the structural analysis input. Following the method established in 
reference 14, the iteration process alternates between loops I and II. The flex- 
ibility of the system allows the user to choose the execution sequencing that 
best fits the problem at hand. A typical mix of the loops is three executions 
of loop I for each execution of loop II. Details of the system organization are 
described in reference 15, and an outline of the major analysis programs involved 
is provided in reference 16. All data transfers showed by the flow chart in 
figure 1 are fully computerized. The iterative process continues until the con- 
vergence criteria for aerodynamic loads, structural mass, and stresses are 
satisfied. 

Structural resizing is executed for each structural component separately 
as if that component, together with the forces acting on it from the neighbor- 
ing components, were extracted from the surrounding structure. In this way, an 
independent problem of structural optimization is formulated with respect to the 
component forces which are treated as invariant external loads. Mathematical 
nonlinear programing is used to solve this problem for each structural component. 
After all components have been optimized, new aeroelastic loads are computed, 
the total structure is reanalyzed, and new component forces are generated for 
use in the next round of individual component optimizations. 

As pointed out in reference 17, where this approach and its positive and 
negative implications are discussed in detail, the disadvantage of such a decom- 
position of the problem is that optimization of components is not tantamount to 
optimization of the whole; or, in other words, minimization of the individual 
component masses does not guarantee minimization of the total mass. This situ- 
ation is caused by the inability to control the load path on the assembled struc- 
ture level; in this regard, the procedure resembles the fully stressed design 
(FSD) algorithm. The advantages lie in the drastically reduced dimensionality 
of each of the individual optimization problems that must be solved and in the 
inherent modularity of the solution algorithm. The consequence of the former 
is that the required computer resources stay within reasonable bounds. The mod- 
ularity permits selection of the resizing algorithm considering generality, effi- 
ciency, and accuracy in isolation from the rest of the problem. Details of the 
component optimization procedure are described in the following section. 
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OPTIMIZATION PROCEDURE FOR SANDWICH PANELS 


The object of the optimization procedure for each individual panel is to 
minimize the panel mass subject to constraints of stress, strain, panel buckling, 
and minimum gage. Panel load consists of the membrane forces Nj^, Ny, and N^y 
shown in figure 2. These forces are held constant for the duration of each cgn- 
secutive optimization procedure for individual panels. For metal panels, the 
design variables are the face-sheet thickness and the total sandwich thickness. 

For composite panels, the variables are the thickness of the composite material 
laid at a given orientation angle (layer thickness) , the orientation angles them- 
selves, and the total sandwich thickness. Construction of a composite sandwich 
panel and the design variables are shown in figure 2. Organization of the pro- 
cedure is depicted in figure 3. The flow chart of figure. 3 corresponds to the 
box marked "RESIZE ELEMBMTS" in figure 1. A panel defined by the design vari- 
ables and subject to the given loads is analyzed (box 1 in fig. 3) for strain, 
stress, and panel buckling. Output from the analysis is used to evaluate con- 
straints and the objective function which is the mass of the panel. Subsequently, 
the results of this evaluation are used by a general purpose optimizer (box 2) 
to calculate a new set of design variables modified to reduce the objective func- 
tion and to move toward satisfaction of the panel constraints. Operations 1 
and 2 are repeated until the constrained minimum mass is reached. Then, the pro- 
cedure is repeated for each next consecutive panel until all panels have been 
optimized. When all panels have new dimensions, the wing structure is reanalyzed 
for aerodynamic loads and internal forces (loops I and II in fig. 1) and another 
round of the individual panel optimization procedure is carried out. 


Evaluation of Constraints 

The optimization algorithm requires that the degree of violation or satis- 
faction of the constraints be stated in the following standard, dimensionless 
form: 


giCYfYa) = 0 


(i e I) (1) 


where gi 


— 1 

ya 


symbolizes the ith constraint, y is a behavior variable 


such as stress, etc., and yg is the corresponding allowable value of y dic- 
tated by the material limits (yield stress, ultimate stress, etc.) or resulting 
from the structural stability considerations, i.e., computation of the critical 
buckling stresses. The purpose of the analysis of the panel being optimized is 
to supply values of y and yg used in equation (1). 


Side constraints .- The side constraint (minimum gage) values for the skins 
are set to provide, as a minimum, fpur plies in each face sheet of the sandwich. 
These plies are associated in pairs with orientation angles (j> and -(|). In 
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most cases (j) is selected to be 45°. Thicknesses associated with other orien- 
tations, that is 0° and 90°, are permitted to vanish but not to become negative; 
therefore, their minimum gage is zero. 

Upper limit on the skin thicknesses is provided indirectly by imposing a 
minimum gage value on the sandwich core thickness as a fraction of total sandwich 
thickness. Without such constraint, the face thicknesses could grow inward to 
the extent of eliminating the core. 

Stress constraints .- Standard, anisotropic membrane analysis (ref. 18) is 
used to produce the stress-strain output for each ply for given edge forces and 
plate geometry. The analysis is simplified by the following assumptions: 

(1) The face-sheet plies are symmetric with respect to the midsurface of 
the sandwich panel. 

(2) The material associated with each orientation is uniformly dispersed 
throughout each face of the sandwich. 

(3) The lay-up is orthotropic (t<jj = t_<j)) although the axes of orthotropy 
need not be aligned with the panel sides (general orthotropy) . As shown in 
figure 2, the orthotropy axes are oriented at an angle y with respect to the 
x-axis coinciding with the panel side. 

The result of assumption 1 is a decoupling of the basic force-deformation equa- 
tion, from reference 18, 


V 


A B 

e 

M 


B D 
L J 

k 


into the equations N = Ae and M = Dk. Assumption 2 permits evaluation of 
the terms of matrices A and D independent of the stacking sequence, and 
assumption 3 simplifies the buckling analysis. 

Solution of the membrane problem represented by N = Ae yields, for each 
fiber orientation, stresses parallel and perpendicular to the fibers a-| and 
02 , the in-plane shear stress ^ 12 # and the corresponding strains e-| , € 2 , 

and £i 2" These values are combined with the allowable limits in constraint 
equations in the standard form shown in equation (1). For example, if a-| ^ 0, 


°1 

1=0 ( 3 ) 

^t°l t 
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where o-|t allowable tension and is an associated factor used to 

modify the allowable stress for selected load cases. A typical example is the 
reduction of the allowable stress level for a 1g cruise condition due to fatigue 
considerations. All numerical information on load factors and allowable stress 
reduction factors actually used in this study is reported in the section on the 
vehicle model. 

In order to prevent matrix failure, the strength constaints are of the form 
given in reference 19; 



where 02g^ = ^ 2 ^ if °2 = ^2a ~ ®2c ^2 ^ This constraint (eq. (4)) 

is a dimensionless inequality equivalent to the standard form represented by equa- 
tion (1). All stress constraints are evaluated for each fiber orientation and for 
each loading case. In addition to these constraints, strain components of a sand- 
wich panel are kept from exceeding specified limits by setting the following con- 
straints: 


^ 1 S 0 (5) 

x,a 


where the subscript x (or y) refers to the panel coordinate that has at least 
one axis aligned with one of the panel sides and the subscript a denotes corre- 
sponding allowable strains. These overall constraints are useful since they pro- 
vide a degree of additional control over the panel strain levels. This control 
is important for at least two reasons: (1) To meet the structure's overall stiff- 

ness requirements, i.e., those dictated by anticipated aeroelastic problems while 
still at the strength design stage and (2) To protect the panel edge members from 
overstress which may occur, for example, in the case of titanium spar caps fram- 
ing a composite panel. (See section "Types of Construction.") 

Buckling constraints .- Buckling of -a plate such as the one shown in f igure 2 
cannot be solved exactly in a closed analytical form because of the general ortho- 
tropy and arbitrary in-plane loads. Since numerical buckling analysis is too 
costly in the repetitive application required in optimization, an approximate 
but fast method is used involving analytical solutions given in reference 20 for 
plates of special orthotropy. The special orthotropy assumption results in 
neglecting the bending- twist coupling terms D^g and D 26 and appears to be 
justified for this type of application by the findings of reference 21 . 
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According to reference 20, the general state of load is separated into the 
cases of biaxial normal stress and shear stress. The measure of buckling criti- 
cality <^ue to biaxial membrane forces is given by the quantity 



where m and n are the number of half-waves of the buckling deformation func- 
tions in the x- and y-dir actions, respectively. 

The plate buckles when Z i ^ , for any combination of m and n. The com- 
bination that minimizes Z is searched for numerically. A conservative approx- 
imation to the critical shear is obtained by the formula (ref. 20) , for an infi- 
nite length strip. 


(Nxy) cr 


[pii (022)^]^'^^ 

b2 


(7) 


where c is given numerically in reference 20 as a function of the parameter 
1 *^ 22 ) 

. The combined action of the biaxial normal forces and the shear is 

°33 

described by the interaction equation, which forms the buckling constraint. 


Z + 


N 


xy 


(Nxy) 


cr 


-1^0 


( 8 ) 


Satisfaction of equation (8) guarantees satisfaction of the following 
special cases; 


Z - 1 S 0 


(Nxy = 0) 


(9a) 


and 


N 


xy 


(Nxy) 


cr 


- 1^0 


(Nx = Ny = 0) (9b) 
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Objective Function 


The structural mass of the panel mp is the objective function. The vari- 
able part of that mass is the sum of face-sheet mass and core mass and is a func- 
tion of the thickness variables. In addition, a constant part is also included 
to represent masses of the bonding and of the nonstructural material (m^jgp and 
mindep) according to the formula 

n 

"p - + PcSItc - tt) + mS + m^epS + min^j^pJ (10) 

i=1 


Optimization Algorithm 

The optimization algorithm, represented by box 2 in figure 3, used for the 
optimization of composite panels is a method of feasible-usable directions, orig- 
inally due to Zoutendijk (ref. 22), coded in a program described in reference 23. 
This algorithm turned out to be efficient, converging relatively independent of 
the number of design variables for the reported application. 


VEHICLE MODEL 

The structural sizing algorithm described in the foregoing section is applied 
to the structural box of the low-aspect-ratio arrow wing of a super sonic- transport 
configuration. The basic details of the finite-element model of that configuration, 
as well as loads, materials, and design variables, are given in the following 
section. 


Finite-Element Model 

The finite-element model is shown in figure 4. It is symmetric with respect 
to the airplane longitudinal center line with one-half consisting of 746 grid 
points and 2141 degrees of freedom connected by 2369 structural finite elements. 
The model represents wing, fuselage, vertical and horizontal tails, wing fins, 
engines, and engine mounts. Nonstructural masses of the engines, leading- and 
trailing-edge devices, and fuel tanks with fuel are included. The wing is 
modeled as follows: . 


(1) Spar and rib caps - rod elements having only extensional stiffness 

(2) Spar and rib webs - quadrilateral panels having only shear stiffness 

(3) Wing covers - quadrilateral and triangular membrane elements, isotropic 
for metal covers and anisotropic for composite covers 

(4) Engine mounts - bar elements having extensional, bending, and torsional 
stiffnesses 
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(5) Wing fins and tail - quadrilateral or triangular plane elements having 
in-plane and bending stiffnesses 

(6) Fuselage frames - rod elements 

(7) Fuselage longerons - rod elements 

(8) Fuselage skin - quadrilateral or triangular membrane elements 

Since the fuselage structure is not included in the resizing process, the finite- 
element model of the fuselage is less detailed than the wing. A box beam is used 
for the fuselage and is composed of quadrilateral membrane elements and shear 
panels preserving the mass distribution of the fuselage and its overall bending 
and torsional stiffnesses. Statistically derived distributed nonoptimum and 
nonstructural masses are accounted for by appropriate coefficients. 


Loading Cases 

From the multitude of loading cases usually considered in airframe design, 
the three cases selected and consistently used in this study were: 2.5g maneuver 

at Mach 1.2, -2g taxi at Mach 2.7, and Ig cruise at Mach 2.7. These cases consti 
tute a minimal set of loadings needed for comparisons of the construction alter na 
tives and for methodology evaluation. The rationale for their selection is as 
follows: (1) the maneuver loads generate the largest bending and shear forces 
near the wing root; (2) negative taxi loads activate buckling constraints in the 
bottom cover of the wing; and (3) constraining the level of stress due to the 
cruise loads is a simplified way to control airframe fatigue life. Because of 
aircraft flexibility and difference in Mach numbers, the maneuver and cruise load 
distributions differ significantly from each other. Therefore, stresses corre- 
sponding to one of these two load cases cannot be obtained by scaling stresses 
caused by the other . 


Types of Construction 

In this study, a structure made entirely of titanium is adopted as a base- 
line design against which the composite construction results are compared. Two 
types of composite construction are considered for the wing: (1) a mixed, or 

hybrid, type in which only the sandwich-panel face sheets are replaced with 
composite material while the remainder of the structure is made of metal and 
(2) a pure composite construction in which all parts, except the aluminum honey- 
comb sandwich core, are made of composite material. 


Material 

The titanium data used are displayed in table I. The allowable stress 
level for cruise is restricted to the value of F^^ to approximately account 
for fatigue requirements. In titanium construction, this value corresponds to 
a notch factor of 4.0. The data for the honeycomb core and core face -sheet 
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bonding are also included in table I. The composite material used in hybrid 
and pure composite construction is graphite-polyimide. 

Graphite-polyimide is a composite material of interest for supersonic- 
cruise applications because of its relatively good retention of mechnical 
properties at elevated temperatures. The material properties anticipated for 
1986 (ref. 2) are displayed in table I. In addition, a reduced allowable strain 
is considered to establish a safe bound on the fiber stress allowables. The 
fiber volume is assumed to be 60 percent throughout. The composite is assumed 
available in a low Young's modulus, high-strength version and in a high Young's 
modulus, low-strength version. In the discussion to follow, use of the low- 
modulus version is assumed unl.ess otherwise indicated. Insofar as the fatigue 
stress allowable for cruise is- concerned, no data similar to those used for 
titanium construction exist for composites. However, for consistency, the same 
ratio of the allowable cruise stress to the design limit stress is used for com- 
posites and for titanium. 


Design Variables and Fixed Parameters 

In the design studies reported in the next section, wing area, planform, 
airplane mass, and structural layout are kept constant. Subject to change are 
the type of construction and the material. Within each of these discrete choices, 
the continuous variables for hybrid and pure composite constructions are selected 
from the set and include sandwich depth ("caliper" thickness) t^,; thicknesses of 
the layers tg, and tgg associated with orientation angles 0°, *|), and 

90°, respectively; and the orientation angles y and 4). (See fig. 2.) Thus, 
the largest number of variables is six since the balanced laminate assumption 
of ^ holds. In order to address the obvious manufacturing requirements, 

grouping of some of these design variables over blocks of adjacent panels is 
also considered. For the case of metal construction (baseline) , the variables 
are sandwich thickness t^ and face thickness t^. 


RESULTS 

Numerical results obtained for the arrow-wing structure are presented in 
this section. The results pertain to three major consecutive phases of the 
study: (1) comparison of construction type and material, (2) optimization with 

the design variables determined most effective, (3) reoptimization with reduced 
number of design variables’ (grouping) to Vssess penalty for uniformity required 
for manufacturing. Results discussed in this report represent a strength design 
in which no attempt is made to control the wing stiffness and flutter. 


Comparison of Construction Type and Material 

Several alternative designs are compared in terms of total sandwich-cover 
structural masses, distributions of the sandwich-cover face-sheet thickness over 
the wing, criticality of the buckling constraints, and stiffness characteristics. 
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Thickness distribution .- The distribution of total thickness t^ is shown 
for the titanium baseline configuration in figure 5(a) and for the composite 
configurations in figures 5(b), (c) , (d) , and (e) . The corresponding masses 
are recorded in table II. The masses represent the total of two wings, includ- 
ing nonoptimum mass. 

In this phase of the study, the design variables are face-sheet thickness 
for the titanium baseline configuration and thicknesses associated with 0°, ±45°, 
and 90° angles of fiber orientation for the composite configurations, with sand- 
wich depth fixed at t^ = 1 in. for all configurations. A high Young's modulus, 
low-strength version of the composite material is used. As seen in figure 5(a), 
the distribution is characterized by the thickness plateaus A and B, which corre- 
spond to the location of the wing fin and engine mount, ^respectively. A similar 
thickness distribution for the hybrid construction can be seen in figure 5(b). 
Maximum thicknesses for the hybrid construction (figs. 5(b) and (d) ) are about 
twice the thickness for titainium construction. However, this is more than com- 
pensated for by the lower density of the composite material, which results in 
the mass savings discussed subsequently. 

Eliminating the titanium caps and changing the shear web material from 
titanium to ±45° composite turns the structure into a pure composite construc- 
tion with the, skin thickness distribution shown in figure 5(c). With titanium 
caps removed, the composite must carry a greater load; therefore, thickness builds 
up (as seen in fig. 5(c)) to values higher by about 25 to 50 percent, . compared 
with the hybrid construction. 

Thickness contours of the individual plies of 90° and ±45° orientation angles 
are shown in figure 5(e) for the LM/H configuration. The contours (fig. 5(e)) 
indicate a characteristic concentration of the ±45° orientation material in the 
area where the load path curves toward the wing root from the direction established 
by the wing sweep angle. Concentrations typical for total skin thickness are 
reflected in the distribution of the 90° orientation (spanwise) material. 

In some instances these concentrations illustrate a well-known instability 
potentially inherent in a fully stressed design approach. Under that approach, 
an overstressed element is stiffened, thus attracting more load which in turn 
requires more stiffness in the next resizing cycle. An example of such thickness 
buildup is designated "region C" in figure 5(d). This numerical instability can 
be checked by imposing upper-limit side constraints on the thickness variables 
in the optimizer. 

Criticality of buckling constraints .- The panels which evolve from the 
optimization as being critical in buckling are shown in figures 6(a) and (b) . 
Comparison with figure 5(d) indicates correspondence of the thickness ridge and 
row A of buckling critical panels in figure 6(a). Removal of the titanium spar 
caps puts all the load on the composite panels, thus increasing very signifi- 
cantly the buckling critical area of a low Young's modulus wing cover. (See 
fig. 6(b) .) The change to high Young's modulus material is effective in reduc- 
ing the number of buckling critical panels to just two panels in figure 6(a) 
and one panel in figure 6(b). 
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structural mass and stiffness characteristics .- Attention is now turned 
to such overall structural characteristics as structural mass and stiffness. 

Total structural mass of two wings includes mass of face sheets, core, bonding, 
ribs, and spars. Both webs and caps are accounted for in the latter, except 
for a pure composite construction in which there are no structural caps. Appro- 
priate nonoptimum mass is added to all structural components. Static stiffness 
can be characterized by wing tip deflection (relative to the wing root) under 
2.5g pull-up maneuver load, and stiffness can be partially measured by the 
natural frequency of the wing first bending mode. 

The three wing characteristics, mass, tip deflection, and fundamental fre- 
quency, are represented by a vertical bar diagram in figure 7 for the B, LM/H, 
HM/H, m/P, HM/P, and LM/H-CA configurations. Comparison of the bars shows a 
significant, 20-percent decrease in structural mass brought about by changing 
from the B configuration to the I2i/H configuration. That benefit is partly off- 
set by an increase in the wing tip deflection and a loss in the frequency of the 
wing first bending mode. 

The loss of stiffness is caused by the fact that the ratio of the allowable 
stress of the composite to that of titanium is higher than the ratio of the 
respective Young's moduli, and is also strongly influenced by changes in thick- 
nesses associated with different orientation angles. Lower strength but higher 
Young's modulus of the HM/H configuration reduces the mass by 13 percent and 
increases stiffness, as manifested by a drop in the wing tip deflection. 

A drastic mass reduction of 56 percent is observed in the LM/P configura- 
tion, but it is countered by an equally sharp increase of the wing tip deflection. 
However , there is a small decrease of frequency because of the attendant decrease 
in mass. These effects are caused by the removal of the titanium caps with their 
contribution to the mass and their limiting influence on the composite strains. 
(See previous discussion on the constraints protecting the titanium caps from 
over stress in a hybrid construction.) 

The HM/p configuration diminishes the mass savings to 49 percent but radi- 
cally improves stiffness, as indicated by significant reduction of the wing tip 
deflection and rise in the wing frequency, the latter being tempered somewhat 
by growth of the structural mass. The influence of the imposition of a conserva- 
tive strain allowable in the LM/H variant (IM/H-CA configuration) is seen to cause 
the mass reduction to drop to 16 percent (from 21 percent noted previously). 

The wing tip deflection decreases slightly from the baseline, but the wing fre- 
quency increases because of reduction in mass. 

In summary, data collected in figure 7 and table II show the maximum mass 
savings to be 56 percent of the baseline mass when the configuration used is the 
low Young's modulus pure composite construction (IM/P) . 


Selection of Design Variables 

The subset of design variables consisting of three thicknesses for orienta- 
tion angles fixed at 0°, ±45°, and 90° has been used so far. Effects of freeing 
the other design variables are now assessed for the Ul/H-CA configuration. 
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Distributions of thicknesses and orientation angles .- Distributions of the 
thicknesses are illustrated in figures 8(a) and (b) , which show contours of total 
thickness of all layers in both faces of the upper wing cover sandwich panels. 
Figure 8(a) corresponds to the layer thicknesses tg/ tgg# sandwich depth 

tg/ and orientation angle <1>. All thicknesses t and t^ are free design vari- 
ables while orientation angle y i-s kept constant at a value of 90°. Figure 8(b) 
corresponds to the set of design variables reduced to tg/ t<j), tgg, and Y 
while freezing t^j and (j) (<1> = 45°); figure 8(c) illustrates the result of free- 
ing the variable t^ and thicknesses tg, t,j), and tgg while freezing the 
angles y and (j) at 90° and 45°, respectively; and finally, figure 8(d) depicts 
the result for all six design variables being simultaneously free. Comparison 
of the results in figure 8 shows significant changes of the cumulative ply thick- 
ness distribution caused primarily by the smoothing effect of the variable t(,. 
This effect is due to the fact that the sandwich depth increase is a very effi- 
cient means of alleviating the panel buckling. 

Individual ply thicknesses and orientation angles (|) are presented in fig- 
ure 9(a) for variable (|> and constant y. Orientation angle is indicated in 
figure 9(a) by arrows for each panel (numerals at each arrow denote the associ- 
ated ply thickness) . The orientation-angle distribution reveals that diagonal 
fibers associated with angle <t> tend to align with the spanwise fibers in the 
rear of the inboard wing box, which is dominated by bending moment. A signifi- 
cant number of angles ^ close to conventional ±45° orientation is found in 
the front upper portion of the inboard wing box, where the load path is turning 
toward the fuselage as it comes from the outboard wing. 

Numerical values of angles y and (j) when both are used as optimization 
variables are shown in figure 9(b). It is apparent that the original setting 
of Y = 0° (°t close to it) is preserved by the optimizer over a great portion 

of the wing area. The only wing areas where it is advantageous to alter y 
in the optimization process are in the vicinity of the wheel well, over the out- 
board wing tip, and near the fin attachment, where the structural layout causes 
abrupt changes of the load path direction. 

Influence on buckling behavior .- The location of the buckling critical 
panels over the wing depends significantly on the choice of design variables, 
as revealed by comparing the different combinations in figures 10(a), (b) , (c) , 
and (d) with figure 6(a). Adding more design variables permits satisfaction 
of the stress-strain constraints with less material and, therefore, makes buck- 
ling significantly more critical; this is illustrated by the greater number of 
buckling critical panels in figures 10(b) and (d) than in figure 6(a). Increase 
in the number of such panels over the lightly loaded front part of the inboard 
wing (fig. 10(c)) indicates that the sandwich depth t^, is decreased to reduce 
the core mass. 

Changes in mass and stiffness characteristics .- Variations of the thickness 
distributions caused by different choices of design variables have a cumulative 
impact on total mass, as indicated by values collected in table III. The table 
is arranged to show progressively lower mass from top to bottom brought about by 
the design variable choices listed. This arrangement reveals a general trend 
of structural mass decrease brought about by greater design freedom, that is, a 
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larger number of design variables which are free to change under the control of 
the optimization algorithm. Thus, a change from "thickness only" to "all six" 
design variables reduces the mass of the sandwich covers (sum of column 3 and 
column 4) by 21 .5 percent and the total mass by 9.4 percent. A comparison of 
mass differences points to the sandwich depth t^, as the most effective vari- 
able in reducing mass. Angle <() consistently causes more change than angle y, 
which is the least effective variable. The effect of 4> is somewhat magnified 
when coupled with t^,, as shown by comparison of relative mass decrements in the 
third and seventh rows of table III. 

The wing-tip-deflection column of table III shows that activating the angle 
design variables reduces stiffness more than structural mass. The variable t^, 
has a similar effect, although it appears to be less detrimental to stiffness 
than the variables <j> and y. The loss of stiffness which accompanies the mass 
reduction can be attributed to the lack of the stiffness constraints in the 
optimization. 

Results presented in this section identify the sandwich depth t,, and orien- 
tation angle (j) as the design variables, besides the ply thicknesses, which 
appear to be most effective under the particular set of loads and constraints. 
However, there are obvious manufacturing constraints that limit the design vari- 
able freedom. It lies within the procedure capability to address, at least par- 
tially, some of these constraints, as described in the following section. 


Variables Constant Over Large Areas of Wing Cover 

In the foregoing discussion, the manufacturing constraints ate set aside, 
temporarily, in order to permit an unrestricted exploration of the benefits that 
may be brought about by freeing those design variables, such as orientation angles, 
which under a conventional approach would probably remain fixed. 

Among the design variable choices considered, the proposition of varying 
orientation angles <1) and y for each panel separately is the least producible 
one. However, inspection of the <t> and y distributions over the wing, such 
as shown in figure 9(b), suggests the possibility of grouping. 

Examples of grouping for variables Y and (j) are illustrated by fig- 
ures 11(a) and 11(b), respectively. Similarly, grouping of sandwich depth t^, 
is depicted in figure 11(c). The effect is that rather large wing areas, e.g., 
almost the whole "glove" region, may be assigned a single value of a particular 
variable. The question of whether or not the groupings shown go far enough to 
render the extra manufacturing cost justifiable can only be answered on a case- 
by-case basis; the point to be made here is simply that such a grouping can be 
made. 


Examples of two grouping alternatives are illustrated in figure 12 for 
both upper and lower wing covers. Grouping 2 (fig. 12(b)) is more restrictive 
than grouping 1 (fig. 12(a)) as it results in a more uniform (more producible) 
distribution of angle <t>- Once a grouping is defined, the optimization is 
repeated with the grouped variables input as constants in the optimization pro- 
cedure. The purpose of this operation is to adjust the remaining free variables 
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to the new situation (new load path, new constraints becoming critical) which 
is created by altering the previous optimization result by grouping. Such 
adjustment by means of reoptimization restores a major part of the mass reduc- 
tion gained by means of the original optimization (ref. 24) . The results of 
the reoptimization are shown in table IV in the form of sandwich cover face- 
sheet total mass for the two groupings portrayed in figures 12(a) and (b) . The 
table lists a 9.3-percent mass reduction gained by freeing 4> in addition to 
thicknesses, when compared with a reference (0°, ±45°, 90°) design in which only 
the thicknesses are free. That mass reduction diminishes to 7.8 percent and 
6.2 percent for groupings 1 and 2, respectively. The results indicate that two- 
thirds of the weight reduction realized by freeing the variable (j> can still be 
recovered when that variable is grouped and frozen. As expected, a more restric- 
tive grouping corresponds to a lower percentage of recovery. 


Summary of Trends 

Results reviewed in the preceding section appear to have design implications 
which can be summarized as follows. 

Structural mass savings of 56 percent can be obtained for the subject airframe 
sized for strength only by changing from baseline titanium to pure composite con- 
struction of a high-strength (low Young’s modulus) material. Mass reductions due 
to hybrid construction and use of composite material having lower strength but 
higher Young’s modulus values are contained in that range. These mass reductions 
correspond to the use of a conventional 0°, ±45°, 90° lay-up and a constant depth 
for the sandwich wing cover panels. 

Further mass reductions on the order of 5 percent are possible if the number 
of design variables is increased to include orientation angles and sandwich depth. 
Most of that increased mass reduction potential is realizable if one addresses 
the producibility constraints by reoptimization with grouping and/or linking. 

These structural mass trends are depicted in figure 13, which shows (from left 
to right) the structural mass of the construction variants discussed in the fore- 
going, relative to the titanium baseline taken as 100 percent. 

Except for the thickness design variables, sandwich depth consistently 
emerges as the variable most effectively influencing structural mass, suggesting 
a design with stepped or tapered-wing sandwich cover depths. The high Young’s 
modulus material option appears to be the most effective means to counter the 
stiffness decrease associated with the change to composite materials. 


Convergence and Computational Requirements 

The method's convergence is illustrated in figure 14 by iteration histories 
for the two configurations B and HM/H. Loops I and II (see fig. 1) are typically 
performed in a ratio varying from 1/1 to 3/1, that is, one to three aerodynamic 
loads updates (loop I) per one structural resizing (loop II) . One complete itera- 
tion is defined as one execution of resizing loop described by the flow chart in 
figure 1 plus the associated load update loop executions. Quantities plotted 
(fig. 14) are structural mass, wing root bending moment, and wing tip deflection. 


16 



The curves indicate a smooth monotonic convergence for both configurations. 

All three quantities, and particularly the root bending moment, stabilize very 
fast in the baseline configuration. Convergence is not as fast for the bending 
moment and tip deflection in the composite construction because of the composite 
wing's greater flexibility, and consequent stronger coupling between the struc- 
tural deflections and aerodynamic loads. The mathematical optimization of each 
panel of the wing cover requires from 8 to 15 objective function and constraint 
evaluations and terminates mostly by the criterion of small consecutive changes 
of the objective function. In summary, the whole iterative procedure appears 
to be well behaved and converges to a reasonable tolerance within five itera- 
tions. Approximately, this consumes as much computer resources as seven repeti- 
tions of a complete aerodynamic load and finite-element structural analyses of 
the subject airframe. 

The procedure depicted in figure 1 is composed of sequential execution of 
several computer programs. The major programs are listed in table V with their 
field lengths and CPU time for one execution for a CONTROL DATA CYBER 175 com- 
puter system. Computer performance indicates acceptable consumption of computer 
resources and generally produces a v;ing optimization within five iterations, CPU 
700-|o sec, and less than 24 hours of turnaround time. 


CONCLUDING REMARKS 

A computerized procedure to generate cross-sectional dimensions and fiber 
orientations for composite airframe structures has been developed. A methodology 
of its application to a wing structural synthesis has been established involving 
the selection of the type of construction, material, and design variables. The 
algorithm unifies computations of aeroelastic loads, stresses, and deflections 
and optimal structural sizing and fiber orientations in an open-ended system of 
integrated computer programs. Aerodynamic loads analysis, finite-element struc- 
tural analysis, and a mathematical-optimization technique are the main components 
of the procedure. Each wing cover panel is optimized separately for minimum mass 
subject to stress, strain, local buckling, and minimum gage constraints.. The 
design variables are ply thicknesses and fiber orientation angles. Experience 
gained in application to a low-aspect-ratio wing of a large supersonic transport 
airframe leads to the following conclusions regarding the procedure: 

1 . A unified iterative method combining structural resizing and aeroelastic 
load computation converges smoothly and monotonically for both metal and composite 
wings. 


2. Use of mathematical programing in an integrated open-ended and modular 
system of computer programs provides flexibility needed for variety of design 
variables, types of material and construction, and constraints; including some 
of the producibility requirements that have to be dealt with in composite struc- 
tures. However, since individual panels are optimized separately, the resulting 
design is expected to be near, but not at, the minimum total mass design. 

3. Computer performance of the method as implemented on the CONTROL DATA 
CYBER 175 computer indicates acceptable consumption of computer resources and 
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generally produces a wing optimization within five iterations, CPU 700-|o sec 
and less than 24 hours of turnaround time. 


In the particular application, structural mass reductions in excess of 
50 percent result when changing from baseline titanium to pure composite con- 
struction. Only about one-third of that mass saving is realized when the change 
is made to a hybrid construction which retains the titanium network of ribs and 
spars. 

Varying the sandwich panel depth, in addition to the layer thicknesses, is 
found to be a significantly more effective means of reducing the structural mass 
than varying the ply orientation angles. Generally, the mass reduction is asso- 
ciated with loss of stiffness. That loss can be reduced by use of a higher 
Young's modulus material. 


Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
January 5, 1979 
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TABLE I.- MATERIAL PROPERTIES 


(Room temperature, 20® C (68° F)] 


Young's modulus, E 

Poisson's 
ratio, \) 

Ultimate tensile 
strength, Ftu 

Ultimate compressive 
strength, F^u 

Stress allowable 
for cruise, Fgc 

Mass density of 
ply material, p 

Fiber 

orientation 

^ 

GN/m2 

ksi 

MN/m2 

ksi 

MN/ra2 



ksi 

MN/m^ 

ksi 

kg/m^ 

Ib/in^ 

Titanium alloy® 

110 

16 X 1Q3 

0;3 

827 

120 

868 

126 

155 

22.5 

4429 

0.16 


Low Young's modulus graphite-polyimide*^ 

138 

7.8 

^bbi 

0:31 

'18 

2033 

115 

295 

16.7 

1998 

113 

WBEm 

351 

51 

m 


Parallel 

Perpendicular 

High Young's modulus graphite-polyimide^ 

276 

12.4 

40 X 10^ 

1.8 X 1q3 

0.29 

.013 

m 

m 

868 

39.3 


176 

25.6 

m 

0.058 

.058 

Parallel 

Perpendicular 


^Minimum gage per face sheet, 0.04 cm (0.015 in.). 

^Properties anticipated for 1986; fiber volume assumed to be 60 percent. 


Mass of honeycomb core, kg/m^ (Ib/in^) 89 (0.0032) 

Mass per unit area for adhesive bonding face sheets to core, kg/m^ (Ib/ft^) . . . 0.036 (0.288) 

Conservative strain allowable for applied load 0.004 

Single ply thickness, cm (in.) 0.0102 (0.004) 
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TABLE II.- MASS^ AND STIFFNESS DATA FOR SEVERAL CONFIGURATIONS 


Configuration 

(b) 


Mass, 

kg ( Ibm) 



Tip 

deflection, 
m (in.) 

Fundamental 

frequency, 

Hz 

Rib and spar 
caps 

Rib and spar 
webs 

Core 

Covers 

Total 

B 

7 814 

4416 

4 755 

10 285 

27 270 

6.5 

1 .182 


(17 228) 

(9736) 

(10 482) 

(22 675) 

(60 120) 

(256) 


LM/H 

7 814 

4399 

4 980 

4 636 

21 829 

o 

• 

00 

.897 


(17 228 

(9698) 

(10 978) 

(10 221) 

(48 125) 

(313) 


HM/H 

7 814 

4347 

4 915 

6 558 

23 634 

3.9 

1 .099 


(17 228) 

(9584) 

(10 836) 

(14 457) 

(52 105) 

(155) 


LM/P 

0 

1733 

4 954 

5 304 

11 992 

10.4 

1 .082 



(3821) 

(10 922) 

(11 694) 

(26 437) 

(409) 


HM/P 

0 

1400 

4 879 

7 531 

13 810 

3.9 

1 .076 



(3087) 

(10 756) 

(16 602) 

(30 445) 

(153) 


LM/H-CA 

7 814 

4382 

4 917 

6 292 

23 406 

6.3 

1 .033 


(17 228) 

(9661) 

(10 840) 

(13 872) 

(51 601) 

(247) 



^Including nonoptimum mass. 

^Hybrid and pure composite configurations are for fixed core depth and ply orientations of 0°, 
±45°, and 90°. 




TABLE III.- FREE DESIGN VARIABLES FOR LM/H CONFIGURATIONa 


Free design 
variables 

Mass, kg (Ibm) 

Wing tip 
deflection, 
m (in.) 

Fundamental 

frequency, 

Hz 

Caps ^ 

Webs 2 

Cote 3 

Covers ^ 

Total 3 

^±45' ^0 

7 814 

4399 

4 980 

4 636 

21 829 

8.0 

0.897 


(17 228) 

(9698) 

(10 978) 

(10 221) 

(48 125) 

(313) 


^0' ^±45' ^0' Y 

7 814 

4410 

4 983 

4 541 

21 748 

8.1 

.896 


(17 228) 

(9722) 

(10 986) 

(10 Oil) 

(47 947) 

(319) 


^ 0 ' ' ^90' 'I’ 

7 814 

4358 

4 993 

4 288 

21 453 

8.3 

.866 

i 

(17 228) 

(9608) 

(11 008) 

(9 453) 

(47 297) 

(327) 


^ 0 ' ^±(t> ' ^90 » Y ' 4> 

7 814 

4374 

5 000 

4 085 

21 274 

• 

CO 

.861 

- 

(17 228) 

(9644) 

(11 024) 

(9 006) 

(46 902) 

(329) 


to f t+45 , tgo / 

7 814 

4418 

3 920 

4 212 

20 363 

8.0 

.897 


(17 228) 

(9739) 

(8 642) 

(9 285) 

(44 894) 

(314) 


tfl f t+ 45 , tgo» 

7 814 

4435 

3 947 

4 095 

20 291 

8.0 

.898 


(17 228) 

(9778) 

(8 702) 

(9 027) 

(44 735) 

(316) 


'-O' '•±(() ' '•90' ^c' 

7 814 

4410 

3 898 

3 831 

19 908 

8.6 

.847 


(17 228) 

(9722) 

(8 594) 

(8 445) 

(43 889) 

(338) 


tO' tgo, 4>i Y 

7 814 

4382 

3 980 

3 607 

19 784 

8.6 

.845 


(17 228) 

(9661) 

(8 774) 

(7 953) 

(43 616) 

(339) 



^All masses include nonoptimum mass. 


ISJ 

u> 




TABLE IV.- RESULTS OF GROUPING DESIGN VARIABLE (j) FOR 
IM/H CONFIGURATION 


Panel grouping 

Cover 

mass 

kg 

Ibm 

No grouping, (j) free 

1396 

3077 

No grouping, (j) fixed at 45° 

1539 

3393 

Grouping 1 (fig. 12(a)) 

1419 

3129 

Grouping 2 (fig. 12(b)) 

1443 

3181 


TABLE V.- COMPUTER RESOURCES 


Program 

Field length 

CPU time, sec 

La) 

SPAR - a finite element 
analysis program 

165 100 

73.7 

Jig-shape computation 

40 100 

1 .3 

Woodward-Car michael 
aerodynamics 

52 100 

30 

Generation of pressure 
coefficients for camber 
and twist 

52 100 

.3 

Airplane-balance 

calculation 

52 100 

.8 

Rib and spar web resizing 

50 100 

.9 

Skin resizing 



120 100 

93.9 


^CONTROL DATA CYBER 175 computer system. 
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Figure 1 Iterative procedures for aeroelastic loads computation (loop I) and 
wing cover resizing (loop II) (from ref. 14). 








Figure 3.- Resizing-procedure flow chart. 
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Buckling critical for low Young's modulus 
Buckling critical for high and low Young's modulus 
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(a) m/H configuration. 


(b) LM/P configuration 


Figure 6.- Buckling critical panels for two configurations. 
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Figure 7.- Mass, tip deflection, and fundamenti 
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Figure 9.- Individual ply thicknesses and orientation angles. 
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Figure 9.- Concluded. 
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Figure 14.- Iteration histories for B configuration and HM/H configuration 
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